A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
2(12.03.3; 12.12.1) 



00 



Is the Universe a fractal? Results from the 
SSRS2 

A. Cappi\ C. Benoist^, L.N. da Costa^-^, S. Maurogordato^ 

^ Osservatorio Astronomico di Bologna, via Zamboni 33, 1-40126, Bologna, Italy 
^ European Southern Observatory, Karl-Scliwarzscliild-Str.2, D-85748 Garching bei Munchen, 
Germany 



0^ : ^ Departamento de Astronomia CNPq/Observatorio Nacional, rua General Jose Cristino 77, 

^ . Rio de Janeiro, R.J. 20921 Brazil 

?H ! * CNRS, UMR 6527; Observatoire de la Cote d'Azur, B4229, Le Mont-Gros, F-06304 Nice 

■ Cedex 4, France 

00 : 

, Received / Accepted 

> ; 

, Abstract. We perform a fractal analysis of the Southern Sky Redshift Survey 2, foUow- 

00 

ing the methods prescribed by Pietronero and collaborators, to check their claim that 
the galaxy distribution is fractal at all scales, and we discuss and explain the reasons 

00 

Q.^ ' of some controversial points, through tests on both galaxy samples and simulations. We 
confirm that the amplitude of the two-point correlation function does not depend on 
^ 'i the sample depth, but increases with luminosity. We find that there is no contradiction 

O ' 

$^ ' between the results of standard and non-standard statistical analysis; moreover, such 
^ ' results are consistent with theoretical predictions derived from standard CDM models 

• • 

J> of galaxy formation, and with the hypothesis of homogeneity at large scale (^ 100 
• I— I 

Mpc). However, for our SSRS2 volume-limited subsamples we show that the first zero- 
?H . 

point of the autocorrelation function ^(s) increases linearly with the sample depth, and 
that its value is comparable to the radius of the maximum sphere which can be com- 
pletely included in the sampled volume; this implies that the true zero-crossing point 
of ^(r) has not been reached. We conclude that the apparent fractal behavior is due to 
a combination of a luminosity-dependent correlation amplitude and the recovering of 
power at larger scales in deeper samples. 
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1. Introduction 

One of the pillars of the standard cosmological models is the large-scale homogeneity 
of the Universe (e.g. Peebles 1993). The standard statistical methods to analyze the 
large-scale structure, as the two-point correlation function ^(r), are described by Peebles 
(1980); they rely on the definition of a mean galaxy density fi, which is meaningful 
only if the assumption of large-scale homogeneity is true. However, at small scales the 
autocorrelation function of the galaxy distribution is positive and can be fitted by a 
power-law, which is a property of a fractal set (see Peebles 1980). On the basis of the 
observational evidence that galaxies are clustered in ever increasing systems, from groups 
and clusters to superclusters, de Vaucouleurs (1970, 1971) presented "the case for a 
hierarchical cosmology"]^. Mandelbrot (1982) developed this concept, suggesting that the 
galaxy distribution in the Universe is fractal, with dimension D = 1. 

In addition, in the last 20 years redshift surveys at increasing depths have revealed 
ever larger structures and voids (see e.g. Davis et al. 1982; de Lapparent et al. 1986, da 
Costa et al. 1994, Vettolani et al. 1997, and references therein). Einasto et al. (1986) found 
evidence that increases with sample volume, but they estimated that it should reach a 
value ^ IQh^^ Mpc for a fair sample of the Universe. Pietronero (1987) stressed that, if 
homogeneity is not reached, the correlation length tq cannot be taken as a measure of the 
clustering amplitude, and suggested a slightly but significantly different definition of the 
autocorrelation function. Adopting this approach, Coleman et al. (1988) reanalyzed the 
CfAl redshift survey (Huchra et al. 1983), concluding that the distribution of galaxies 
was fractal to at least 20 h^^ Mpc. On the other hand, Davis et al. (1988) found that tq 
increased as Tq'^, and not linearly with the depth as predicted for a simple fractal (see also 
Maurogordato et al. 1992). Others advocated the need for a multifractal approach (e.g. 
Balian & Schaeffer 1989; Martinez & Jones 1990; Martinez et al. 1990); this is however 
another issue and we will not discuss it: see for example the review by Borgani (1995) 
and references therein. 

The first redshift surveys sampled relatively small volumes, and the reality and nature 
of correlation amplitude variations with depth remained an open -and much debated- 
question. Could these variations simply reflect fluctuations due to local structures? Were 
they a consequence of the fractal distribution of galaxies? Or were they an indirect 
consequence of the dependence of clustering on galaxy luminosity? 

^ The idea of a hierarchical Universe has a long history, which dates back to the XVIIl"' 
century (see E. Harrison, 1987, Darkness at night, A Riddle of the Universe, Harvard University 
Press, Cambridge); in this century, Fournier d'Albe and Charlier were the first to propose a 
hierarchical (now we would say fractal) model of the Universe (see Mandelbrot 1982). 
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It should be pointed out that even at a scale 1000 Mpc we do not expect a 
perfect homogeneity, as COBE has found evidence of anisotropy in the CMB radiation 
at a level ST/T ~ 10~^ (Smoot et al. 1992), a value necessary and sufficient, at least in 
some standard cosmological models, to explain the formation of the observed structures. 
The existence of very large structures in the Universe implies that the galaxy or cluster 
autocorrelation function must be positive to a scale corresponding to the size of these 
structures; but this is not necessarily inconsistent with the measured value of the galaxy 
autocorrelation length ro ^ 5 — 6 Mpc for L galaxies (as claimed by Pietronero 
and collaborators), as larger structures have also a lower contrast, and the correlation of 
luminous matter may be significantly amplified relatively to the underlying dark matter 
correlation function (Kaiser 1984; Bardeen et al. 1986). Therefore, we expect that the 
galaxy and cluster distribution will not be perfectly homogeneous even at large scales. 
The claim for a fractal Universe is obviously much more stronger than that: it implies 
that there is no convergence to homogeneity and that it is not possible to define a mean 
density n of the Universe. 

In the last years Pietronero and his collaborators applied statistical indicators which 
do not imply a universal mean density on an ever increasing number of catalogs (see 
e.g. Di Nella et al. 1996, Sylos Labini, Montuori & Pietronero 1996, Montuori et al. 
1997), claiming evidence for a galaxy fractal distribution at all scales. Their results are 
impressive, as recently stressed by Coles (1998), but they appear to be at variance with 
the results of other groups who analyzed the same catalogs with standard indicators. It 
is clear that the situation is unsatisfactory, despite many articles and even public debates 
on the subject (see Pietronero et al. 1997; Davis 1997; Guzzo 1997); while the strongest 
support to large-scale homogeneity comes indirectly from the high level of isotropy of the 
CMB radiation and from two-dimensional catalogs (Peebles 1996), there is still confusion 
on the interpretation of the available three-dimensional data, mainly due to the difference 
in the statistical indicators, which do not allow a direct and quantitative comparison of 
the results. 

Therefore we analyzed the Southern Sky Redshift Survey 2 (SSRS2; da Costa et al. 
1994) following the approach of Pietronero and collaborators, in order to independently 
check their claims and to answer to their criticisms (Sylos Labini ct al. 1997, hereafter 
SMP) about the work of Benoist et al. (1996, hereafter Paper I). In section 2 we discuss 
the apparent dependence of ro on sample depth; in section 3 we describe the different 
statistics, and their relations; in section 4 we present and discuss our fractal analysis of 
the SSRS2 in comparison with the standard analysis; in section 5 we show that theoretical 
predictions derived from the standard CDM model can reproduce our results, and are 
therefore consistent with a homogeneous Universe; our conclusions are in section 6. 
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2. Does ro depend on Scimple depth? 



If the Universe is homogeneous at large scales, then for any fair sample the statistical 
properties of the galaxy distribution should be the same; in particular, the autocorrelation 
function of any given class of galaxies should not depend on the sample depth. On the 
contrary, if the Universe is a simple fractal at all scales, the mean density in a sphere of 
radius r centered on a randomly chosen galaxy will vary as 

n(r) oc (1) 

where D is the so-called correlation dimension. 
The autocorrelation function will vary as 

^(r) oc r(3-^) - 1 (2) 
and the correlation length ro will increase with the sample depth: 

ro = [(3-7)/6]i/^i?. (3) 
where 7 = 3 — Z) and Rs is a measure of the sample size. 

In practice, there is the problem to define a sample not dominated by a few prominent 
structures, and to discriminate between a dependence of ro on the sample depth or on the 
absolute luminosities of the galaxies, as deeper volume-limited samples extracted from a 
magnitude-limited catalog will progressively include more luminous galaxies. The most 
recent results appear to confirm the original claims by Hamilton (1988) and Borner et 
al. (1989) about a real dependence of ro on the luminosity. In particular, in Paper I we 
analyzed the SSRS2 and showed that ro has a significant increase for L > L^, galaxies. 
SMP criticize these results and point out that a value of ro as large as 16 Mpc, as 
found in our Paper I for the M < —21 galaxies, has never been found in previous surveys 
as the CfAl. We stress however that a) the previous, shallower surveys do not contain 
enough galaxies with M < —21 galaxies which, being rare, require a larger volume; b) as 
shown in figure 5 of Paper I, the results obtained by Hamilton on the CfAl arc consistent 
with the new SSRS2 results -up to M = —20.5, representing the maximum luminosity 
reachable in a volume-limited sample large enough for a statistical analysis in the CfAl-. 
Notice that the CfAl and the SSRS2 arc limited respectively to to = 14.5 and to = 15.5, 
so volume-limited samples containing galaxies with the same absolute luminosity have 
different depths; their consistency at a fixed absolute magnitude is another argument in 
favor of luminosity segregation and against a fractal distribution. However, SMP criticize 
our Paper I also because we tested the independence of the correlation length from the 
sample depth only for a few subsamples; moreover, they notice that the results obtained 
by Sylos Labini & Pietronero (1996) on the Perseus-Pisces catalog show a systematic 
dependence of ro on the sample depth, and not on luminosity. Indeed the same test 
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Table 1. Volume-limited subsamples from the SSRS2 
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Table 2. Volume-limited subsamples from the SSRS2 with a limiting apparent magni- 
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performed by SMP on the Perseus-Pisces is possible also for all the subsamples of the 
SSRS2. In figure (|l|) we present ro with its bootstrap error as a function of the absolute 
magnitude for volume-limited samples of the SSRS2 selected with apparent limiting 
magnitudes of 14.5 and 15.5 (see table 1 and 2). Notice that the distance limits for the 
volume-limited samples are calculated assuming iJp = 100 km s^^ Mpc^^ and taking 
into account the K-correction, and they are conservatively fixed to include the galaxies 
of all the morphological types (for more details see Paper I). 

As discussed in Paper I, the less deep subsamples (including galaxies with a faint 
absolute magnitude) may be affected by the dominance of local structures. From figure 
^, there is indeed marginal evidence for sampling variations in subsamples including 
galaxies with absolute magnitude fainter than Af = —19, while for deeper subsamples 
galaxies with the same absolute magnitude have the same value of tq within the errors, 
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Fig. 1. Correlation length tq for different volume-limited subsamples extracted from 
the SSRS2, assuming an apparent limiting magnitude mum = 14.5 (open triangles) and 
miim = 15.5 (solid triangles) 

independently of the sample depth. This consistency internal to the SSRS2, in addition to 
the external consistency between subsamples extracted from the CfAl and the SSRS2, is a 
further and robust confirmation that the correlation amplitude depends on the luminosity 
and not on the sample depth. This result is obviously in disagreement with the claims by 
SMP, based on subsamples of the Perseus-Pisces supercluster, which show indeed very 
large differences between the correlation length ro of galaxies with the same absolute 
magnitude but in volumes at different depths (see their figure 40) . SMP assume a constant 
error Aro = 0.5 Mpc; while a more realistic error would have been preferable, it 
cannot be a large underestimate given the size of their subsamples. 

While there is no a priori reason to prefer one catalog to another, the type of in- 
formation which can be extracted depends on the particular selection criteria of each 
catalog. In this case, contrarily to the "anonymous" CfA or SSRS redshift surveys, the 
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name Perseus-Pisces redshift survey reminds us that it is a region of the sky selected 
because of the presence, at about 5500 km/ s, of a large filamentary structure (Giovanelli 
& Haynes 1991). While many interesting statistical results can -and could- be obtained, 
their interpretation must take into account this peculiarity of the survey. For example, 
Guzzo (1997) discusses the effect of the Perseus-Pisces supercluster on the galaxy counts, 
questioning the interpretation of Pietronero et al. (1997). As concerning the variation of 
To with sample depth, we notice that a volume-limited sample including all galaxies 
brighter than -19 has a depth of ~ 5000 km/s at an apparent limiting magnitude of 
m = 14.5, and a depth of ~ 8000 km/s at an apparent limiting magnitude of m = 15.5. 
In the case of the Perseus-Pisces redshift survey, this means that in the first case we 
include the void in front of the supercluster but not the supercluster itself, while in the 
second case we include the whole supercluster. This might at least partly explain the 
discrepancy with our results; therefore a fractal galaxy distribution does not appear as 
the most plausible explanation for the systematic variations of ro in the Perseus-Pisces 
redshift survey. 



3. Statistical tools 

As we mentioned in the introduction, the statistical indicators used to assess the fractality 
of the Universe are different from -but not independent of- the standard ones. In this 
section we will clarify the relations between the two sets of indicators. 

A simple dimensionless measure of galaxy clustering is the two-point correlation 
function ^(r), which is defined by the joint probability 5P of finding a galaxy in a 
volume element dVi and another galaxy in a volume element dV2 at separation r: SP = 
n^[l + ^{r)]dVidV2, where we take into account only the magnitude of the separation 
but not its direction, assuming as usual that the galaxy distribution is a stationary 
point-process (see Peebles 1980). The two-point correlation function ^(r) of a sample is 
measured generating a catalog of randomly distributed points in the same volume, with 
the same selection function, then counting the number Ngg{r) of distinct galaxy-galaxy 
pairs at separation r, and the number Ngr{r) of galaxy-random pairs at separation r 
(Davis & Peebles 1983): 

where nc is the mean galaxy density and ur the density of points in the random 
catalog. Hamilton (1994) suggested a slightly different estimator of ^(r) less dependent 
on the density, and best suited to measure ^(r) at large separations, which we adopted 
in Paper I, and which we use also in this paper. This is what we will call in the following 
the "direct" estimate of ^(r). 
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The galaxy two-point correlation function is well approximated by a power-law, 
^(r) = (r/ro)^'', where 7 ^ 1.8, at least at scales less than ^ IQh^^ Mpc: this implies 
that the correlation dimension £' = 3 — 7 is ^1.2. 

As discussed in the previous section, if the Universe is a fractal, it is impossible to 
find a "universal" value of tq , because larger samples will have larger voids and a larger 
r-Q. For this reason Pietronero (1987) has defined the conditional average density T{r): 

2 ^ I i-r+dr 

rM-^E^/ nir.+r)dr (5) 

which measures the average density at distance r from an occupied point. r(r) is in 
principle related to ^(r): 

C(r) - T{r)/n ~ 1 (6) 

One implication (trivial but important) of relation (^ is that, if ^(r) is a power- 
law, r(r) is not, and vice-versa. Moreover, there is a practical difference in the way of 
estimating the two functions (see e.g. Coleman & Pietronero 1992). T{r) is calculated only 
within the sample limits, while ^(r), as apparent from equation (^, is computed using 
random catalogs to correct for volume boundaries, thus implicitly assuming homogeneity. 
Even if the Universe is not a fractal, this can be a problem when the sample size is still 
small relatively to the scale where homogeneity is achieved; as a consequence, an indirect 
estimate of r(r) obtained from ^(r) through equation (^ may be wrong. 

^(r) is however a simple statistical measure of clustering; other useful functions rep- 
resenting variously defined integrals of ^(r) can be directly measured on a sample. 

For example, the mean number of objects Np(r) in a volume V centered on a randomly 
chosen object, is 

Np{r)= [1 + ^{x)]dV ^ An x^[l+^{x)]dx^TiV + 4:nJ3 (7) 
Jo Jo 

where Jsir) is the so-called J3 integral: 

Mr) = / x^a^)dx (8) 
Jo 

AnJ^ represents the mean number of excess galaxies around a galaxy within a radius 

r. 

Another important integral quantity is the volume-averaged two-point correlation 
function ^(r): 

^r) = ^ j^^{r,r2)d\,d\2 (9) 

Contrarily to ^(r), ^(r) is usually computed through counts in randomly placed spher- 
ical cells, taking into account only the spheres totally included within the sample bound- 
aries. No random catalog is used, analogously to the estimate of r(r). 
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The equivalent of Jsir) and ^(r) in the fractal approach is the integrated conditional 
density r*(r) defined as: 

o nr 

T*{r) = - — - / ATrx'^Tix)dx (10) 







which "gives the behavior of the average density of a sphere of radius r centered 
around an occupied point averaged over all occupied points" (Coleman & Pietronero 
1992). 

We notice that r*(r) is related to Jsir): using relations (H), (^) and (10) we have: 



r*{r)^n[l + ^Mr)] (11) 
which represents the integral counterpart of equation (^. 

As in the case of ^(r), J3{r) is also computed using random samples, i.e. assuming 
homogeneity beyond the sample limits, so equation (|l^) cannot be used to obtain T* (r) 
from J3. 

Assuming now that f (r) is a perfect power law, we have the following relations: 

^3(.) ^ 3(r^e(o (12) 

and 

^"^^^ ^ (l-7/3)(l'7/4)(l-7/6)^('^) ^''^ 
We have then: 

r*(r) = n[l + 2^(1 - 7/4)(l - 7/6)^»] = n[l + (J^^l (14) 



Equations (|l2[-|l4|) are useful because they allow us to estimate ^ (r) from the integral 
indicators; in the next section we will use them, with equation (^), to compare the 
different functions we have defined above. 



4. Fractal analysis of the SSRS2 

A preliminary test -and the easiest one- in order to compute the fractal dimension of 
the galaxy distribution can be performed counting the number N(r) of objects within 
a distance r from us: from equation (^, we expect N{r) oc . The results for the 4 
volume-limited SSRS2 samples with absolute magnitude limits M < — 18, M < —19, 
M < —20, AI < —21, are shown in figure ^. A formal least-squares fit to the data gives 
D = 2.11 ± 0.27, D = 2.59 ± 0.21, D = 3.34 ± 0.18, D ^ 3.40 ± 0.35; the errors were 
computed through bootstrap resampling of the original samples. These results apparently 
tell us that from the shallowest to the deepest sample we have a change in the fractal 
dimension, with convergence to homogeneity. Even considering the uncertainty associated 
to these values, we can conservatively state that such results are more consistent with 
homogeneity at scales ~ 100 Mpc, than with a fractal Universe. 
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Fig. 2. Integrated number of galaxies as a function of conioving distance r for 4 volume- 
limited subsamples of the SSRS2 {M < -18, M < -19,A/ < -20,M < -21). The solid 
lines are least-squares fits to the data vi^ith the law N{< r) oc , with respectively 
D = 2.11 ± 0.27, D = 2.59 ± 0.21, D = 3.34 ± 0.18, D = 3.40 ± 0.35. 



We have also computed r(s) and r*(s) for the same 4 volume-limited samples. In 
what follows we will use redshift space estimates. First of all, we are interested here in 
a relative comparison between different indicators, and not in determining a true value; 
second, we want to compare our results with those of Pictroncro and collaborators, always 
performed in redshift space; finally, the SSRS2 is not affected by large peculiar motions 
(see Willmer et al. 1998). 

The results are shown in figure ^ A least-square fit in the small-scale range (less 
than 10 h^^ Mpc) gives a slope -f = 3 — D 1, i.e. D 2. At larger scales, the slope 
fluctuates from sample to sample, and within the same sample, but does not show a clear 
plateau (as should be expected if homogeneity was reached); moreover, positive values 
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Fig. 3. r*(s) (panel a) and T{s) (panel b) for the SSRS2 volume-limited samples in- 
cluding galaxies with M < —18 (triangles), M < —19 (squares), M < —20 (hexagons), 
M < —21 (circles). The solid line has a slope 7 = — 1 (i.e. D = 2). 

of r(s) and T*{s) arc measured to larger scales for deeper samples. These results are 
consistent with those found by Pietroncro and collaborators on other catalogs. 

Does all this mean that the Universe is fractal, and that the hypothesis of large-scale 
homogeneity is falsified? 

First of all, we have to take into account the behavior of N{r) we have discussed 
above, which is consistent with homogeneity at large scales. r(s) and r*(s) are more 
general tests, as they average on all the galaxies, while N{r) is centered on our position; 
nevertheless, with N{r) we can test larger scales (> 100 Mpc), while r(s) and r*(s) 
are limited to ~ 40 Mpc. At this scale we still have a positive correlation function 
(see below) and we should not expect to measure D = 3; therefore what we can infer 
from all these tests is that there is still power at a scale of ~ 40 Mpc, but no evidence 
for significant power at ~ 100 Mpc. 
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A further problem for the fractal model is the discrepancy between the slope of ^ (s) , 
which implies D ^ 1 (for the SSRS2 see our results in table 1 of Paper I), and that of 
r(s) and r*(s), implying D ^ 2. SMP claim that people fit ^(s) in the "wrong range" 
around ro, where ^(s) would deviate from a power-law. We notice that people usually 
fit ^(s) where it can be reasonably approximated by a power-law, excluding both small 
and large scales; therefore the explanation of SMP (see their figure 9) does not seem 
plausible. 

In interpreting these results, it should perhaps be taken into account an effect which 
is trivial but which we point out for sake of clarity. In figure ^ we show ^(r) and V{r) 
assuming that ^(s) is a perfect power-law with 7 = —2, and points are obtained assuming 
an error of 20%. It is obvious that the quantity V{r) cx [1 + i{r)] can be approximated 
by 2 power-laws; in this case, in the range 2-10 Mpc it has a slope 7 ~ —1.3 or 
D — 1.7. At larger scales, where ^(r) becomes very small, F becomes constant. However, 
in real samples ^(s) deviates from a power-law at small and large scales (see below), and 
the regime where T{r) is constant may not be reached. This simple exercise tells us that if 
^(s) is a genuine power-law, there is a range around tq where T{s) can be approximated 
by a power-law with a flatter slope, which is what is observed. 

We will check now the internal consistency between standard and non standard statis- 
tics, examining the behavior of the functions described in the previous section. In figure 
^, for the volume -limited subsamplc including all galaxies with M < — 20, we show the 
values of ^(s) estimated directly using the Hamilton (1994) formula, from T{s) through 
equation (j^), and from the integral functions Jj,{s), ^(s) and r*(s) using equations (p^), 
(p^), ( p^ (i.e., assuming that ^(s) is a perfect power law). In the standard scenario, 
all these estimates applied on the same sample should in principle measure the same 
quantity, and the only differences should be due to the biases introduced by each method 
(this is why no error bars are shown in the figure). 

Figure || shows that there are two main effects: at small scales, there is a systematic 
difference between the estimates of ^(s) obtained from the integrated functions and the 
two estimates obtained directly or from F(s); at large scales, there is a systematic dif- 
ference between the estimates based on random catalogs and the other ones which only 
take into account the volumes fully included within the sample limits. 

The direct estimate of ^(s) is consistent with that computed from F(s) to about 
lO/i^^ Mpc. At small scales peculiar motions affect ^(s), and the number of pairs at 
small separations becomes very small; therefore we expect significant deviations of ^(s) 
from a power-law, and this explains why the indirect estimates of ^(s) from r*(s), ^(s), 
and Jz{s) based on the assumption of a perfect power-law, differ systematically from the 
other two estimates. 
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Fig. 4. ^(r) (dashed line) and r(r) = n[l + ^(r)] (solid line) assuming ^(r) = (ro/r)^, 
with ro = 5h~^ Mpcand n= 1. Open and filled triangles simulate a measure of the above 
two functions assuming a 20% error. 

At scales larger than ^ lOh^^ Mpc the estimates of ^(s) obtained from r(s), T*{s) 
and ^(s) decline more rapidly than the direct estimate of ^(s) and that obtained from 
J3(s). As the former indicators take into account only the volumes fully included within 
the sample limits, they are limited to smaller scales -their maximum scale corresponds 
to the radius of the largest sphere which can bo included in the sample-; moreover, at 
large scales the randomly placed spheres have a significant overlap, and are no more 
independent, thus giving a biased (lower) estimate of the variance. It can be noticed that 
the estimate of ^(s) from T{s) drops somewhat more rapidly than the estimates from 
r*(s) and ^(s), because the latter are derived from integral measures, thus are smoother 
and have a better signal to noise ratio. It is also clear that the estimate of ^(s) from J^is), 
which is an integrated measure and is based on random catalogs, can reach the largest 
scales, and appears to be consistent with a power-law extending to ~ 50 Mpc. The 
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Fig. 5. ^(s) derived with different statistical tools for the M = —20 subsample of the 
SSRS2. Open triangles: direct standard estimate; filled triangles: computed from F; open 
squares: computed from J3; open hexagons: computed from ^; crosses: computed from 

direct standard estimate of ^(s) is noisier and oscillates; it drops beyond 10 h^^ Mpc, in 
a way similar to the integral estimates, then has a second oscillation, following the J3(s) 
estimate. 

We can conclude from the above discussion that there is no disagreement between the 
results obtained with standard and non-standard indicators, when taking into account 
the range of validity of each one of them. 

In figure 1^, for the volume-limited subsamples from A/ = —16 to M = —21 (step 0.5) 
we show the correlation length tq, the first zero-point i?^ of ^(s) (i.e. the scale where ^(s) 
becomes negative), and the radius of the largest sphere totally included in the sample, 
Rmax, as a function of the sample depth Dg. The values of tq with their bootstrap error 
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Fig. 6. The correlation length tq (filled triangles), the zero-point Rc (open triangles), 
and the maximum radius Rmax (open squares) as a function of the sample depth D, for 
volume-limited samples of the SSRS2 (from M < —16 to M < —21, with a step of 0.5 
magnitudes). 

bars are the same as in figure |^, while the values of Rc simply correspond to the last bin 
where ^(s) is positive, and their error bars indicate the bin width. 

It is trivial that Rmax oc D^. The correlation length 7'o is not proportional to the 
depth, contrarily to what expected for a fractal; the less deep subsamples are probably 
not "fair" , and ro may be underestimated (even if a luminosity effect for low luminosity 
galaxies cannot be excluded) ; then ro is more or less constant in an intermediate range of 
luminosities, and increases rapidly (more rapidly than the depth) for samples including 
galaxies brighter than Af». This is the luminosity effect we have shown in Paper I (see 
also figure |l|). The most striking result of figure ^ is that the value of the zero-point R^ 
clearly increases linearly with the depth, and that its value corresponds in a striking way 

to Rmax- 
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This quantitative correspondence between Rc and Rmax means that the standard 
estimate of the correlation function ^(s) is in fact reliable only at scales smaller than the 
largest sphere contained in the sample: it is at this scale which ^(s) is forced to become 
negative, then starts to oscillate. We stress that this effect should be taken into account; 
often an oscillating correlation function is interpreted as evidence of a real characteristic 
scale, but if this scale corresponds to the Rmax of the sample taken into account, the 
unavoidable conclusion (as in the case of our subsamples) is that the real zero-point of 
^(s) has not been reached. 

5. Comparison with theoretical predictions 

Are the results of our fractal analysis consistent with standard models of galaxy forma- 
tion? 

The behavior of F and T* can be checked by using N-body simulations, which are 
based on standard Friedmann-Lemaitre cosmological models, and assume homogeneity 
at large scales. We used the Adaptive P^AI code of Couchman (1991) to evolve a standard 
Cold Dark Matter universe. We started from the initial linear power spectrum of Bond & 
Efstathiou (1984) with Qh = 0.5. The box size was 100 Mpc, with 64^ particles. We 
analyzed the simulation at a time when its correlation function was roughly comparable 
to the observed one, diluting it to approach the observed galaxy densities. Figure ^ shows 
^(r) and F*(r) for 2 subsamples extracted from the simulation: the first with a box size 
L = 25 h-^ Mpc, the second with L = 50 h-^ Mpc. The best-fit slope for F* is 7 = -1.0, 
i.e. D = 2 in the first case; 7 = —0.8, i.e. D ^ 2.2, in the second case, ^(r) shows the 
typical cutoff when approaching the radius of the largest sphere which can be contained 
in the sample: in fact, the deeper sample allows to recover power at larger scales. These 
results are comparable to our results on the SSRS2, and also to those presented by 
SMP. We have already mentioned that a problem for a simple fractal Universe is the 
discrepancy between the fractal correlation dimension implied by ^(s), I? = 3 — 7 1, 
and that derived from F(s) or F*(s), D ^ 2. We have also noted that such a discrepancy 
cannot be due simply to an improper fitting of ^(s) (as claimed by SMP) in a range 
where it deviates from a power-law. 

Now we have shown that the observed slopes of ^(s) and F(s) can be reproduced 
by a standard CDM N-body simulation: this confirms the implications of our example 
in figure ^, and the consistency of both ^(s) and F(s) with gravitational clustering in a 
Universe homogeneous at large-scale. 

We have not addressed, of course, the question of the variations of tq. In our simula- 
tion, To at any fixed time has a unique, constant value; when analyzing deeper subsamples 
of our simulation, tq rapidly converges to its true value (no significant difference in fact 
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Fig. 7. r*(r) and ^(r) measures for one realization of a CDM universe. Filled and open 
triangles represent respectively ^ and F* measured in a box with L = 25 h^^ Mpc; filled 
and open squares represent respectively £, and F* measured in a box L — 50 h^^ Mpc. 
The solid line has a slope 7 ~ —0.8; the dashed line has a slope 7 = —2.0. 

is found in the value ro measured in a box of 50 h^^ Mpc and in another at 75 h^^ 
Mpc), while the scale of the zero-point increases significantly. 

As we have previously shown, in the real Universe the correlation length rg of galaxies 
does depend on their luminosity (or, in other terms, there is evidence of luminosity 
segregation). Our simple answer to the criticism by SMP -"the authors [...] have never 
presented any quantitative argument that explains the shift of with sample size^^~, is 
that shallow samples are affected by the presence of local structures, as the discrepancy 
between the Perseus-Pisces and the other surveys reminds us; moreover, when a sample 
is not deep enough, ^(s) can be artificially truncated, thus affecting the estimate of rg; 
finally, and it is our main quantitative argument, from figure H it is clear that does 
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not depend on sample depth, but on galaxy luminosity (as shown in paper I), for L > 
galaxies. 

On the other hand, there are various reasons why we should expect luminosity segre- 
gation. It has been shown that the correlation function can be amplified as a consequence 
of a "statistical bias" (Kaiser 1984), which can explain in a natural way the cluster auto- 
correlation function. As concerning galaxies. White et al. (1987) used N-body simulations 
to show that more massive galaxies are more clustered; they called this phenomenon "nat- 
ural bias". Recent semi-analytical models (Mo & White 1996, Kauffmann et al. 1997) 
predict this bias if galaxies form in sufficiently massive halos, with M > M*; in this 
case, the trend is similar to that observed in the SSRS2 (see figure 5 of Paper I), even 
if they do not succeed in fitting both the low and high luminosity samples. This is not 
surprising, given all the theoretical uncertainties, particularly on the dark halo-galaxy 
relation. 

In addition to the luminosity -dependent correlation amplitude, there arc other 2 
important points which should be considered when analyzing the galaxy distribution: the 
K-correction and the cosmological model (sec Guzzo 1997; Scaramella et al. 1998). While 
in our analysis of the SSRS2 samples we have fixed an — 100, qo — 0.5 model and 
systematically applied K-corrections, Pietronero and collaborators usually perform their 
analysis in Euclidean space and without K-correction. As concerning the cosmological 
model, for a fractal Universe one should find and apply the appropriate formulae for 
an isotropic but inhomogeneous Universe (unless one wants to question the validity of 
General Relativity). On the other hand, the K-correction is an empirical correction to 
the magnitude of a galaxy, as in a fixed wavelength range we observe different parts 
of the spectrum at different redshifts. The effect is still small in the deepest sample of 
the SSRS2; if we use Euclidean distances and do not apply the K-correction, from the 
relation A^(< r) oc we find D = 3.28 instead of £* = 3.40. It becomes critical in deeper 
surveys as the ESO Slice Project (Vettolani et al. 1997; Zucca et al. 1997), as shown by 
Scaramella et al. (1998), who do not confirm the results obtained by SMP on that sample 
as well as on Abell clusters. 

Moreover, the K correction cannot be ignored when interpreting deep counts; the 
claim by Sylos Labini et al. (1996) that deep comits are consistent with a fractal dis- 
tribution can be reversed, concluding that the agreement with the observed count slope 
would be lost if at least the K correction had been applied (not to speak of cosmologi- 
cal and evolutionary corrections: they are model-dependent, but we know that galaxies 
evolve). 

Why do the choice of an Euclidean space and the absence of a K-correction give 
results more consistent with a fractal distribution? There are essentially two main rea- 



A. Cappi et al.: Is the Universe a fractal? 19 

sons: first of all, Euclidean distances become systematically larger than comoving ones 
at increasing redshift: this means also that one will have larger voids and structures. Sec- 
ond, when applying the K-correction we can recover galaxies that otherwise would not 
have been included in a volume-limited subsample; if the K-correction is not applied, a 
systematically larger fraction of galaxies will be lost at larger distances. The combination 
of the two effects lowers the density at increasing redshift, and flattens the slope of the 
\og{N) — log(r) relation, thus following the same trend of a fractal distribution. 

6. Conclusions 

In this paper, we have carefully examined the claim that the Universe is a fractal, analyz- 
ing the SSRS2 and using the same statistical approach of Pietronero and collaborators. 
Here are our main results: 

— the SSRS2 has a luminosity dependent correlation amplitude; 

— results obtained with the standard and fractal approach give consistent results; 

— concerning large-scale homogeneity, observational evidence is -for the moment- con- 
sistent with theoretical predictions derived from standard CDM models of galaxy 
formation, and does not require a fractal Universe; 

— on the other hand, the first zero-point of the correlation function we measure in 
SSRS2 volume-limited subsamples increases linearly with the sample depth, and has 
approximately the same value of the radius of the largest sphere which can be included 
in the sample. This implies that the zero-point of ^(s) is beyond 40 Mpc, a lower 
limit set by very luminous galaxies (see Benoist et al. 1996; Cappi et al. 1998) and 
consistent with results from clusters (e.g. Cappi & Maurogordato 1992); on this point, 
we agree with SMP that measurements of ^(s) beyond a scale corresponding to the 
radius of the largest sphere included in the sample are not reliable. 

We conclude that the dependence of the galaxy correlation function on the galaxy 
luminosity and the power still present at large scales (> 40 Mpc) can satisfacto- 
rily explain those effects interpreted by SMP as evidence of an inhomogeneous, fractal 
Universe. 
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